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Abstract 

The decay process a — » tttt at high temperature and density and its relation with chiral phase 
transition are discussed in the framework of the NJL model. The decay rate for the process a tttt 
is calculated m the whole temperature and density region. The contribution of the final state pion 
statistics to the decay rate is discussed. The maximum decay rate at different chemical potential is 
computed. Finally, we investigate the relation between the starting point of the decay process and the 
critical point of the first-order chiral phase transition. 
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1 Introduction 

It's generally believed ||l| that there are two kinds of QCD phase transitions in relativistic heavy ion 
collisions. One of them is related to the deconfinement process of quarks, and the other is about chiral 
symmetry restoration. When we investigate chiral properties, a meson has to be included While the 
(T mass and its width are very large in the vacuum, around the chiral phase transition point a becomes 
light and its width is also very small Therefore there will be numerous cr's in the early stage of a 
relativistic heavy ion collision if the system undergoes a chiral phase transition. When the temperature 
and density of the system drop down, a mass increases but pion mass keeps as a constant. When the 
mass condition rria- — 2x71^^ is reached, the process a ~> tttt starts immediately. From these considerations, 
if chiral phase transition happens in the system, the numerous cr's around the phase-transition point will 
affect remarkably the final state pions in magnitude and in momentum distribution |Q| . 

Nambu-Jona-Lasinio(NJL) model 1^ is often used to investigate spontaneous chiral symmetry break- 
ing in the vacuum and chiral symmetry restoration at finite temperature and density Although the 
model has no confinement mechanism in the vacuum and therefore can not be used to study deconfine- 
ment 1^ , it describes chiral properties very well . In this article we calculate in the framework of the 
NJL model the decay rate for a — s- tttt and investigate its relation with the chiral phase transition. 



2 The decay rate for a ^ tttt 

The Lagrangian density of SU(2) NJL model is 

Lnjl = ^(»7^9^ - mo)^ + GH^/j^P)^ + {^i-f5T^b)% (1) 

where ^ and ip are the quark fields, r is the isospin generator, G is the coupling constant with dimension 
GeV~'^, and toq is the current quark mass. 

Taking quark propagator in the mean field approximation and meson propagator in the RPA approx- 
imation ||], we get the lowest order Feynman diagram for the decay process cr — > tttt sketched in Fig. |l|. 
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The decay rate for the process at finite temperature and density is 



8tt ml ^ ■ 2 

where the triangle factor ^cttttt is defined by |^ 

A_(r,^) = 4m,N^NfJj^ — X 

8(g • P)^ - (2m2 + Aml)q- p + /2 - 2mlE^ 
{ml - 4£;2)((m2 - 2q- p)^ - mlE^) ' 

In Eq. T and fi are temperature and bayonic chemical potential respectively, m^ and are a mass 
and TT mass, t/o- and are the coupling constants of qqa and qqi: interactions. The temperature and 
density dependence of these quantities can be found in Ref. [|. feix) = (e^/^ - 1)-^ in Eq. (|) is the 
Bose-Einstein distribution function for the final state pions. In Eq. (^, A^c = 3 and Nf = 2 are the quark 
color and flavor degrees of freedom, frix) = (e^/-^ + is the Fermi-Dirac distribution function for 



(3) 



quarks, Eg — ym'^ + q^ is the quark energy. The temperature and density dependence of quark mass 
mq can also be found in Ref. [||. 




Figure 1: Feynman diagram for the process a — > tttt to the lowest order. The solid lines denote quarks 
and anti-quarks, the dashed lines denote mesons. 



3 The relation with the chiral phase transition 

We first discuss the medium effect for the final state pions in the decay process. The temperature 
dependence (/i = 0) of the process cr ^ tttt was considered in Ref. The difference between Ref. Q 
and Eq. (H) is in the consideration of the statistics for the final state pions. It can be proved that 
the temperature and density effect of two final state mesons can be represented as a statistical factor 
(1 + /i + /2). In the case of cr —> tttt, the correct statistical factor is therefore (1 + 2/b(^)) in Eq. (^, 
but not (1 + in Ref. ||^. This difference in the final state statistics will result in an essential 

change in the behavior of the decay rate at the critical point of chiral phase transition. 

The parameters of the NJL model in chiral limit can be determined by the quark condensate density 
< uu >=< dd >— {~0.25GeV)^ and the pion decay constant f-^ = 0.093GeV in the vacuum(T = 0,/i = 
0), they are the coupling constant G = 5.02GeV~^ and the quark momentum cutoff A = 0.653Gey. 
The decay rate To-tttv as a function of temperature T (/j, = 0) is displayed in Fig. ^. The decay begins 
and reaches suddenly the maximum decay rate at the critical temperature Tc = IQOMeV. The solid 
line represents the decay rate with statistical factor (1 + 2/b) and the dashed line represents that with 
(1 + /b)^. It can be found in Fig. that, there is almost no difference between the two factors at low 
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Figure 2: The decay rate with different statistical factors in the chiral limitfii = 0). 
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Figure 3: The decay rate in the whole T-fi region. 



temperature. However, when the temperature approaches to the critical point of the second-order chiral 
phase transition, the first statistical factor makes the decay rate finite at the critical point, but the later 
one makes it infinite. From this we can conclude that the chiral phase transition leads to a strong pion 
enhancement, but the singularity in the decay rate comes from the wrong consideration of the medium 
effect for the final state pions. By further computation, we can get the analytic expression of the decay 
rate at the critical point, 



FiT-TrTr (-^ ^ ) 



SttTc 
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(4) 



The pion mass is not zero in the real physical world. The explicit chiral symmetry breaking, namely 
mo 0, must be considered. In this case, the three parameters in the NJL model can be determined 
by the quark condensate, the pion coupling constant and the pion mass — O.lSAGeV in the vacuum, 
they are G = 4.936*6^"^ A = 0.6535Gey and mo = 0.005Gey. Fig. | shows the decay rate in the whole 
temperature and density region beyond chiral limit. In the vacuum, the decay rate is about 0.09GeV. 
Keeping T=0 and increasing /i, the decay rate behaves as a constant until fi is close to the critical point 
of the first-order chiral phase transition. The decay rate decreases slightly a.t n — Q.3GeV, then increases 
and reaches the maximum value at the critical point. Finally it jumps down to zero at the critical point. 
In the high density region where the first-order phase transition exists, this jump retains at the chiral 
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Figure 4: The maximum decay rate at different chemical potential fj,. 
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Figure 5: The solid line denotes the first-order chiral phase transition line around endpoint E, the dashed 
line denotes the threshold point for the decay at different fi. 

phase transition points. Whether the decay rate jumps to zero will be discussed in a while. When 
the chemical potential fj, is not so high that first-order phase transition can't happen(/i < 0.327GeV), 
with increasing temperature the decay rate rises gradually from its vacuum value to the maximum at 
some temperature, and this temperature decreases with increasing chemical potential. After reaching 
the maximum, the decay rate drops down rapidly but continuously to zero at the threshold temperature 
defined by (T, n) = 2m7r (T, ^J■) ■ The fact that the maximum value of F is not in the vacuum but near 
the threshold point is attributed to the contribution of the final state pion statistics which is important 
only at high temperature. 

For each chemical potential /x, we can extract a maximum decay rate Tmax from r(r, fi). Fig. |4| shows 
the Tjnax as a function of /i. When chemical potential is very low or very high, namely in the limit of high 
temperature or high density, Tmax is larger than the values in the middle region. However, the variation 
is not big, we can consider T^ax as a constant approximately. 

In the case of explicit chiral symmetry breaking, the second-order chiral phase transition at high 
temperature in the chiral limit becomes a smooth crossover, while the first-order phase transition at 
high density still exists with a new critical point In Fig. |^, the solid line represents the critical 

temperature Tc of the first-order chiral phase transition at different chemical potential /it, the point E is 
the endpoint of the phase transition. As discussed above, the threshold point of the decay a tttt is 
defined by the mass equation m^{Td, Hd) = 2m7r(Td, ^d)- The question is if the threshold point coincides 
with the critical point, or if Td equals to Tc- 

The dashed line in Fig. || represents the threshold temperature Td at different /i. It is easy to see 
that, Td is larger than Tc in the density region near the endpoint E. The difference between Td and Tc is 
maximum at E, and decreases with increasing chemical potential. At some chemical potential it vanishes. 
Therefore, in the first-order phase transition region near the endpoint E, the decay rate jumps at the 





4 



phase transition point, but it does not jump to zero but to a finite value. After this jump the decay rate 
drops continuously to zero at Td- In the rest first-order phase transition region where the critical point 
and the threshold point coincide, the decay rate jumps to zero. Although the decay a tttt happens 
before the phase transition around the endpoint E, our numerical calculations show that only about 1% 
of (t's decay before the phase transition. 



4 Conclusions 

Whether the chiral phase transition happens in relativistic heavy ion collisions can be manifested by the 
number and momentum spectra of final state pions. In mean-field approximation, the chiral symmetry 
restoration changes mainly the masses of particles , and in turn the thresholds of particle reactions. Since 
the different temperature and density dependence of a and tt, the process a — > tttt at finite temperature 
and density changes the properties of final state pions remarkably. Around the chiral phase transition 
point a becomes very light and is easy to be generated. These produced cr's will decay into two pions 
when the decay condition m^r — 2m7r is satisfied. When we consider the medium effect for the final state 
pions correctly, the enhancement will not become unphysically infinite. The almost chemical potential 
independent maximum decay rate means that the chiral properties reflected by a — > tttt can be measured 
in relativistic heavy ion collisions at different energy. While in the chiral limit the threshold point of 
a —>■ TTTT coincides with the critical point of the chiral phase transition, in the case of explicit chiral 
symmetry breaking, the decay may happen before the chiral phase transition. Our numerical results 
showed in this paper depend certainly on the NJL model, but we believe that the above qualitative 
conclusions are model independent. 
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